On the phase of magneto-oscillations in graphite 
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The problem of Dirac fermions in graphite subject to a perpendicular magnetic field is studied. 
We show analytically that the weak inter-layer interaction between the graphene sheets leads to 
anomalies in the Shubnikov-de Haas and de Haas- van Alphen magneto-oscillations governed by the 
orbits around extremal cross-sections of the graphite Fermi surface. The calculation of the Landau 
plot performed within a four band continuum model reveals that magneto-oscillations are aperiodic, 
except of the case of vanishing inter-layer interaction at the H point of the graphite Brillouin zone. 
Also for all other orbits along the H-K-H edge the magneto-oscillations are only asymptotically 
periodic in the quasi-classical limit, with the phase corresponding to massive fermions. 
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Graphite is a layered material composed of weakly cou- 
pled two-dimensional (2D) graphene sheets formed by 
hexagonally arranged carbon atoms. In 2004, a single 
sheet of graphene was prepared from three-dimensional 
(3D) graphite by micro-mechanical cleavage The dis- 
covery immediately attracted attention of the solid-state 
physical community, as the electrons in graphene obey a 
linear energy dependence on the wave-vector and be- 
have like massless relativistic particles - Dirac fermions 
(DFs). In the seminal papers [H, 0, [1], the Shubnikov- 
de Haas (SdH) magneto-oscillations in graphene were 
found periodic in an inverse magnetic field, simi- 
larly as in the case of a 2D gas of massive Schrodinger 
fermions (SFs), but with the phase shifted by tt. The 
shift, which was clearly demonstrated by the Landau plot 
of magneto-resistance oscillations, is due to the existence 
of the zero-energy Landau level (LL), shared by electrons 
and holes. For the same reason, the anomalous quantum 
Hall effect with a half-integer instead of integer quanti- 
zation was observed in mechanically ex-foliated samples 
• This is considered as the most direct evidence of DFs 
in graphene. In 2006, important technological progress 
was achieved. The epitaxial graphene was grown on the 
single-crystal silicon carbide by vacuum graphitization 

The discovery of DFs in graphene has resulted in re- 
newed interest in bulk graphite. In a series of papers 
[1, 0, HI, the spectral analysis of SdH and de Haas- van 
Alphen (dHvA) oscillations was employed to determine 
the phases of two series of magneto-oscillations observed 
in graphite. In papers 0, [13] an attempt was made to 
relate the phases to the topological Berry phase, which is 
acquired by fermions moving around close orbits. Based 
on their analysis the authors of Refs. [f| 0, @ came to a 
conclusion that one of two groups of oscillating carriers 
corresponds to DFs. Recently a paper [ll[ was devoted 
to a careful, mostly experimental investigation of SdH ef- 
fect in graphite, and doubts about the observation of DFs 
using magneto-transport measurements were expressed. 

Here the problem is treated from the theoretical point 
of view. We construct a Landau plot for the model 



Hamiltonian developed by Slonczewski, Weiss and Mac- 
Clure (SWM) [H, [[!, Oj and compare the result with 
Landau plots for SFs and DFs, described below. 

The energy spectra of 2D SFs and DFs in a zero mag- 
netic field can be written as 
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E D (k) = ±hv F k, 
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where k 



Ik 2 , + k 2 , m* is the effective mass of SFs and 

Vp is the Fermi velocity of DFs. The positive and nega- 
tive branches of the conical DFs spectrum correspond to 
the electrons and holes, respectively. 

In a magnetic field, the spectra of SFs and DFs are 
quantized into the LLs as follows: 



E, 



En = ±j2h\e\v 2 F Bn, (2) 



where uj c — \e\B/m* is the cyclotron frequency and the 
index n = 0, 1, 2, • • • . In the case of SFs the equidistant 
LLs lie above E = for any finite B, whereas in the DFs 
case the lowest electron LL is shared with the highest 
hole LL located exactly at E = 0. 

The Eqs. © are consistent with the Onsager-Lifshitz 
quasi-classical quantization rule 



A Q { E F )= 2 -^.(n + 1 Q), Q 



= S,D, 



(3) 



where A®(Ep) = nkp is the area of the SF or DF 
Fermi circle, calculated with the Fermi energy Ep and 
the Fermi wave- vector k F taken from Eqs. |T]). We get 



s _ 



1/2 for SFs and ~/ D = for DFs. 



T 

Magneto-oscillations observed in SdH and dHvA ef- 
fects are controlled by oscillations of the density of states 
(DOS). It is well known that the DOS on the Fermi level, 
g{Ep), can be expressed as an imaginary part of the re- 
solvent G{z) = (z- H)- 1 , 
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lmTrG(E F + i0). 



(4) 
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For simple diagonal Hamiltonians of SFs and DFs given 
by Eqs. ([!]), we get 



2ttH 



(5) 



It follows from Eq. §5§ that the DOS reaches maxima 
at magnetic fields B n for which the LLs cross the Fermi 
energy Ep. A Landau plot, i.e., the plot of the inverse 
magnetic fields 1/B n versus the level index n is a stan- 
dard tool used to determine the frequency and phase of 
magneto-oscillations. For SFs and DFs we arrive to 
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where By 
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n + 7 Q , 



and Bfi> 



(6) 



E F /(2h\e\v F ) are 



the oscillation frequencies, in agreement with the quasi- 
classical expression obtained from Eq. ([3]) , 



B® = hA Q /2n\ 



(7) 



It is clear that the positions of maxima of SF and DF 
oscillations differ by a half of the period, i.e., by n in 
terms of a phase factor. 

In graphite, the inter-layer interaction of Bernald- 
stacked graphenes adds a fc z -dependence to the electron 
energy spectrum and a 3D Fermi surface (FS) is formed 
close to the H-K-H edge of the hexagonal Brillouin zone 
(BZ). As mentioned above, the graphite 3D electronic 
structure is described by the semi-empirical SWM Hamil- 
tonian, which employs seven nearest-neighbor tight- 
binding (TB) parameters 70, 71, ... 75, A, and the value 
of the Fermi energy, Ep. Previously, the model param- 
eters were fitted to various optical and transport exper- 
iments 15j. Recently, their values are continuously re- 
fined by fitting to the experimental data [ll| and/or to 
the results of first-principles numerical simulations of the 
graphite band structure [Til ]. 

Among the seven SWM parameters, the parameter 73, 
which controls the trigonal warping of the FS, brings a 
numerical complications in the case of nonzero magnetic 
field. When 73 is taken into account, the magnetic-field- 
SWM Hamiltonian has an infinite order and must be di- 
agonalized numerically [l7j . Fortunately, its influence is 
not too strong for k z far from the H point of the BZ 
and energies close to Ep [ijj • To facilitate our analytical 
treatment, we prefer to use a simplified Hamiltonian H 
with 73 neglected, in the form introduced in the Mac- 
Clure's paper (l4| . 

The choice of 73 = yields isotropic equienergetic con- 
tours. The FS of graphite consists of elongated elec- 
tron and hole pockets located near the points K and 
H, with /c z -dependent circular cross-sections. In the 
quasi-classical limit, two extremal cross-sections define 
two series of magneto-oscillations and two quasi-classical 
frequencies Bq and Bq for electrons and holes, respec- 
tively. Analytical solutions for the FS cross-sections can 
be found, e.g., in Ref. [l5j . 
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FIG. 1: (color online) a) Contours of the graphite Fermi sur- 
face, b) a dimensionless parameter 5 as a function of fc z , c) 
the parameters B± as functions of k z . 



Within these approximations an expression similar to 
Eq. ([6|) can be obtained for any fc z -dependent cross- 
section of the 3D graphite FS. To do so, we need to find 
the poles of the resolvent G(z) — (Ep — H + iO)^ 1 . In 
other words, we should solve the secular equation de- 
rived from the simplified Hamiltonian H of Ref. [3] for 
B n , i.e., we should find the roots B n of the secular poly- 
nomial. 

The solution yields a formula for the inverse magnetic 
fields as a function of the level index, n, in a shape 



B? 



± 



n(n + l)5(k z 



B n 



(8) 



where the three coefficients B±(k z ) and 5(k z ) can be de- 
rived from the SWM model and the value oi Ep. The 
^-dependence originates from cos (k z c/2) which appears 
in the inter-layer TB SWM parameters, c/2 denotes the 
inter-layer distance in graphite. Obviously, BjL(k z ) / B ni 
as given by Eq. ((5J), are not periodic in 1/B. 

An expression for the Landau plot corresponding to 
Eq. ijSJ) can be written in a form 



Br, 



'it 



(k z ). 



(9) 



In this equation 7^ n (k z ), defined by the right-hand side 
of Eq. ([5]), is no longer a constant describing the oscilla- 
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FIG. 2: (color online) 7±,„(fc z ) for maximum cross-sections of 
the electron and hole pockets approximated by y/n(n + l)—n. 
The inset shows results of the linear approximation of first N 
terms (n = 1, 2, • • • , N) of \/n(n + 1). The relative deviations 
from the quasi-classical frequency and phase are shown. 



tion phase, but a variable which depends on the LL index 
n. 

The dependence on n is most pronounced for high mag- 
netic fields, i.e., for small n. In the low-magnetic-field 
limit, with a large number of LLs below Ef, we can write 
n(n + l)5(k z ) > 1/4 and yjn(n + 1) -> n + 1/2. Then 
7±„(fc z ) — > 1/2, and we can conclude that the charge 
carriers in graphite behave, at least as far as the phase is 
concerned, as the SFs. 

Only when we can completely neglect the inter-layer 
interaction, as at the H point of the 3D BZ, k z = ir/c, 
cos (fc 2 c/2) = 0, A = and S(k z ) — > 0, we get 



I ± i, 

2 2 



(10) 



a result which corresponds to DFs. 

Here the Landau plots of two series of magneto- 
oscillations, which can be observed in graphite, are 
considered as most interesting. We constructed them 
based on the parameters of the SWM model taken from 
Ref. [HI]. The explicit expressions for B± and S(k z ), 
which appear in Eq. (jSJ), will not be presented here. In- 
stead, their fc z -dependences are shown in Fig.[TJ together 
with the contours of the FS. The positive sign applies 
in the formula © for k z from the electron region of FS 
and the parameter is equal to the quasi-classical fre- 
quency Bq. Similarly, the negative sign should be taken 
for k z from the hole pockets, where B_ becomes equal 
to Bq. For both electron and hole extremal orbits the 
parameter 5 is large and, consequently, the right-hand 
side of Eq. ((5]) can be approximated by \Jn(n + 1). The 
marked difference is only for n = 0. The Fig. [2] shows to 
what extent the accuracy of this approximation is rea- 
sonable. 

The expression for B±/B n is not periodic and the 



question arises how many LLs must be resolved to reach 
the linear dependence on n, i.e., the quasi-classical limit 
n+1/2. The result of the linear approximation of a model 
curve yjn(n + 1) is presented in the inset of Fig. [2], where 
the relative deviations from the quasi-classical frequency 
and phase are shown. It turns out that both frequency 
and phase are underestimated if we took into account 
only limited number of n, 1 < n < N. Less oscillations 
are necessary to get close to the quasi-classical frequency 
than to obtain a reasonable approximation for the phase. 
This may explain the differences found between the ex- 
perimentally determined phases of samples with different 
mobilities, which are determined from different number 
of oscillations resolved. 

Note that the above Landau plot can be derived if we 
approximate the energy spectra of electrons and holes by 
the formulae 



El 



Hu e c y/n(n + l), E h n = hu*y/n(n + l), (11) 



where uj e c and uj^ are the quasi-classical cyclotron fre- 
quencies corresponding to extremal electron and hole or- 
bits. 

Two series of magneto-oscillations discussed above cor- 
respond to d 3> 1 and to the maximum cross-sections of 
the electron and hole pockets. There are another two 
extremal cross-sections, till now not reliably resolved in 
the transport experiments, located around the H point of 
the BZ where 6 < 1. In spite of the fact that, according 
to Ref. [l7j . the parameter 73 has a qualitative influence 
on the LLs structure near this point, it is at least of the 
theoretical interest to study the behavior of the 7± n (k z ) 
for model with 73 neglected. 

For a given k z the dependence of the energy bands 
on k is hyperbolic in a zero magnetic field. While 
near the maximum cross-sections there are broad min- 
ima/maxima of bands which are similar to parabolas for 
small k, near the H point the hyperbolas are very sharp 
and with a shape close to the Dirac cone. Therefore, 
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FIG. 3: (color online)The dependence of 7_,„(fc z ) on the LL 
index n for k z close to the H point of the BZ. 



4 



a smooth transition of 7± ra (fc z ) from SFs to DFs is ex- 
pected. Fig. [3] reveals that the behavior is more compli- 
cated. 

In the SWM model the fc-dependence of the zero-field 
energy at the H point (k z = 0) is given by 



E G = ^±J(^J +Vvlk\ (12) 

which is not equal to the Dirac cone for a finite A. Never- 
theless, 6 = implies that 7_ n (k z ) is a constant equal to 
0, as for the DFs. This is in agreement with the Landau 
plot constructed from the energy spectra in a magnetic 
field, which according to Ref. [14j have a simple analytic 
form 



E^-±^-j + 2h\e\v F Bn. (13) 

On the other hand, the parameter 8 equals 1 for Ep 
crossing the Ei band of the SWM model. According to 
Eq. © this leads to J? >n (k z ) = 1/2, as for the SFs, in 
spite of the fc-dependence not so close to parabolic one 
as the for the maximum cross-sections, 5 3> 1, where the 
energy spectra can be approximated by Eqs. (fTTj) . 

The field dependence of the corresponding LLs ranges 
from that described by Eq. (fl"3|) at the H point, which is 
close to \f~B~ characteristic for DFs, to the linear depen- 
dence on B typical for SFs for extremal electron and hole 
orbits, as given by Eqs. (fTTj) . 

It follows from the above discussion that the hyperbolic 
fc-dependence of the zero-field electron energy bands, 
which changes considerably depending on the value of 



k Zl yields aperiodic magneto-oscillations when a mag- 
netic field is applied. An exception are two k z in the 
neighborhood of the H point of the BZ. We assume that 
this conclusion is at least qualitatively correct, as one can 
hardly believe that this is just the neglected 73 which 
yields the magneto-oscillation aperiodicity. 

There is another potential reason for deviations from 
the magneto-oscillation periodicity. Unlike the optical 
experiments which involve electrons with energies below 
and above the Fermi energy, the SdH and dHvA magneto- 
oscillations reflect only the properties of electrons with an 
energy equal to Ep. Our treatment is based on the as- 
sumption that the Ep is a constant. This is not quite 
correct as the carrier concentration is a constant and not 
Ep, which should oscillate as a function of B. This can 
be important for lowest LLs in high mobility samples and 
was considered as a single source of oscillation aperiod- 
icity in Ref. [ll|. 

In conclusion, we have found that the magneto- 
oscillations in graphite are only asymptotically periodic 
in the quasi-classical limit, with the phase correspond- 
ing to massive fermions. The quasi-classical limit can 
be reached only exceptionally for samples with very high 
mobility and at very low magnetic field. Therefore, the 
determination of the oscillation phase in samples with a 
limited number of resolved LLs below Ep is not a reliable 
tool for distinguishing between DFs and SFs in graphite, 
due to the aperiodicity of the magneto-oscillations in a 
standard quantum regime. 
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